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Lecture  6 


CELLULAR  FLAMES 

We  shall  now  exanine  the  left  stability  boundary  that  was  uncovered 
in  lecture  5  in  our  discussion  of  ifEFs  (figure  5-3).  The  boundary  is 
associated  with  instabilities  leading  to  cellular  flames,  i.e.  flames 
whose  surfaces  are  broken  u?  into  distinct  luminous  regions  (cells)  separated 
by  dark  lines.  Each  line  is  a  ridge  of  high  curvature,  convex  towards  the 
burnt  gas.  For  a  nominally  flat  flame  these  cells  are  very  unsteady, 
growing  and  subdividing  in  a  chaotic  fashion;  but  curvature,  for  example, 
can  make  them  stationary. 

The  most  striking  manifestation  of  cellular  instability  is  the 
polyhedral  flame,  into  which  the  conical  flame  on  a  Bunsen  burner  can 
suddenly  transform.  The  conical  surface  splits  into  triangular  cells 
forming  a  polyhedron;  for  a  five-sided  flame  the  appearance,  from  above, 
is  much  like  that  of  the  Chrysler  emblem  (figure  1).  The  dark  wedges 
between  the  white  triangular  cells  correspond  to  sharp  ridges;  the  dark 
central  region  corresponds  to  a  tip  with  strong  curvature.  Figure  2  gives 
a  sketch  of  a  five-sided  flame,  derived  from  a  photograph  in  Smith  & 

Pickering  (1929)-  Polyhedral  flames  are  often  stationary,  but  can  spin 
rapidly  about  the  vertical,  making  several  revolutions  per  second. 

We  shall  discuss  chaotic  and  stationary  cellular  flames,  including 
polyhedral  flames,  in  the  framework  of  the  weakly-nonlinear  theory 
pioneered  by  Sivashinsky.  The  constant-density  approximation  will  be 
used  throughout ,  although  perturbations  of  it  will  be  admitted  in  two 
places . 

1.  Chaotic  Cellular  Structure. 

The  nonlinearity  associated  with  the  left  stability  boundary  will 
be  weakest  in  the  neighborhood  of 

-6.1- 


-6.2- 

I  =  -1 ,  k  =  0  , 

a  possible  bifurcation  point;  accordingly  ve  focus  our  attention  there  by 
” — King 

I  +  1  =  0(e),  k  =  0(/e),  o=  0(E2),  (2) 


where  e  is  a  small  positive  parameter  that  will  be  found  to  represent  the 
amplitude  of  the  disturbance.  The  relative  ordering  of  2.+1  and  k  is 
suggested  by  the  parabolic  shape  of  the  stability  boundary,  while  the  order 
of  a  follows  from  the  limiting  form  (5*^9)  of  the  dispersion  relation 
(5.39)  as  e  0.  This  determines  the  growth  rate  of  the  most  important 
Fourier  components  (the  unstable  ones)  of  the  disturbance  when  fc+1 
is  small.  In  terms  of  any  scalar  F  that  represents  the  disturbance  field, 
the  dispersion  relation  is  equivalent  to 


Ft 


-(1+H)F  =  0. 

yyyy  yy 


(3) 


For  1+1  <  0,  this  equation  predicts  unbounded  growth.  Bifurcation 
(with  weakly  nonlinear  description)  is  possible  if  nonlinear  effects,  not 
yet  taken  into  account,  limit  this  growth.  We  shall  first  give  a  heuristic 
argument  to  determine  these  effects  and  then  substantiate  the  result  by 
formal  analysis.  The  argument  consists  in  recognizing  that  equation  (3) 
is  actually  a  formula  for  the  wave  speed,  and  modifying  it  appropriately. 

In  this  connection,  suppose  that  F  determines  the  location  of  the  flame 

5T.2?t  3.5 


x  =  -t  +  eF; 


(U) 


sceea 


sheet  is 


1+eW  +...  with  =  -F^ 


(5) 
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ar.d  equation  (3)  becomes 


W,  =  U?  -  (l+i)F  . 

i  yyyy  yy 


(6) 


Tnis  formula,  determines  the  deviations  of  the  flame  speed  from  its  adiabatic 
value  of  1  due  to  the  reaction-diffusion  effects  that  are  triggered  by  the 
distortion  of  the  flame  front. 

Now,  equation  (5b)  is  mere  kinematics,  valid  in  a  linear  theory  only'. 

The  exact  relation  between  flame  speed  and  displacement  is 


1-eF, 


W  = 


A  O 

/l+e  F 


=  1-eF  -?e2F2+... 
t  y 


(7) 


and,  for  disturbances  with  wave-numbers  of  the  magnitude  (2b),  the  nonlinear 
2  2 

terms  je  F  is  comparable  to  the  linear  term  eF^.  This  suggests  that  the 


nonlinear  generalization 


w 


1 


(8) 


should  be  used  in  the  formula  (6)  and,  when  e  is  purged  from -the  result¬ 
ing  equation  by  writing 


1+1  =  -e ,  n 


/ey,  t 


e2t 


(9) 


(in  accordance  with  the  ordering  (2)),  we  find 


o 

F  +  ?F  +  !+F  +  F  =  0. 

t  n  nnnn  nn 


do) 


Note  that  this  equation  hold  for  i  <  -1. 

Substantiation  of  this  result  requires  a  systematic  asymptotic 
development  in  which  x  is  replaced  by  the  coordinate 


n  =  x  +  t  -  eF(n,t) 


(11) 


in  the  governing  equations  (h ,2h ,25) ;  thus, 

3/3x  =  3/3 .n,  3/3y  =  -e3^?  3/3n+e?3/3n,  3/3t  =  (l-e3F  )3/3n+e23/3t  (12) 

n  T 

when  y,  t  are  replaced  by  n,t.  The  normal  derivative,  required  for 
the  jump  conditions  (1*. 27-29),  is 

(l+5e3F2)3/3n-e3F  3/3n  (l^) 

n  n 

to  sufficient  accuracy.  Perturbation  expansions  in  £  are  now  introduced 
for  T,  h,  and  F,  leading  to  a  sequence  of  linear  problems  for  the  T-  and  h- 
coefficients  as  functions  of  n,n,  and  t.  These  are  to  be  solved  under  the 
requirements:  T^  =  =  T^  =. .  .=  0  for  n  >  0;  conditions  as  n  -*■  -®  are 

undisturbed;  and  exponential  growth  as  n  -<-+<”  is  disallowed.  The  problems 
are  overdetermined,  but  only  at  the  fourth  (for  T^/n^)  a  solvability 
condition  required,  namely  equation  (10)  for  the  leading  term  in  F. 

For  two-dimensional  disturbances  of  the  flame  sheet,  the  basic  equation 
is 

Ft+  s(VF)2+  kVkF  +  V2F  =  0.  (1U) 

Discussion  for  both  one-  and  two-dimensional  disturbances  has  been  limited 
no  numerical  computations.  The  solutions  obtained  display  chaotic  variations 
in  a  cellular  structure,  resembling  the  behavior  of  actual  flames.  Figure  3 
clearly  shows  the  ridges  that  separate  the  individual  cells. 

2.  Ffrect  of  Curvature. 

Equation  (12 }  is  a  balance  of  small  terms;  it  may  be  modified  to  account 
for  any  additional  physical  process  whose  effect  is  also  small.  Hydrodynamic 
effects  can  be  incorporated,  for  example,  if  the  density  change  across  the 
flame  is  appropriately  small  (because  of  small  heat  release),  and  this 
provides  important  insight  into  the  role  of  Darrieus-Landau  instability 
in  actual  flames.  Equation  (qo  )  is  replaced  by 
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p 

F  +  +  1*?  +  F  + 

t  e  "  nnnn  nn 


“  F  (n,t)  _ 

-D= -  dn  =  0 

n-n 


with  y 


o-l 


2ne 


3/2  ’ 


(15) 


tie  requirement  on  the  heat  release  being  y  =  0(l):  the  expansion  ratio 

3/2 

a  can  only  differ  from  1  by  0(e  ).  Numerical  integration  shows  that  the 

new  (integral)  term  is  destabilizing;  an  even  finer  structure  is  super¬ 
imposed  on  the  chaotic  cellular  pattern  obtained  without  it.  We  shall 
not  consider  the  topic  further,  since  Sivashinsky  (1983)  has  recently 
discussed  it  in  detail.  Instead,  we  shall  examine  a  much  simpler  effect, 
that  of  curvature. 

Consider  a  line  source  of  mixture  supporting  a  stationary  cylindrical 
flame;  the  radial  speed  of  the  efflux  is  taken  to  be 


u  =  R/r, 


(16) 


where  R  is  an  assignable  constant.  The  corresponding  solution  of  equations 
' - . 2-  ,25 )  ,  the  jump  conditions  (-.27-29),  and  the  boundary  conditions  at 
zr.s  origin  and  infinity  show  that  the  flame  is  located  at  r  =  R  and  that 


',+Y,(r/a) 


R 


(tYfrRln(R/r)/RH_1 


h  = 


(17) 


for  r  <  R. 


Note  that  the  speed  (lo)  is  1  at  the  flame  location,  an  unexpected  result. 

Apparently  the  effect  of  flame  curvature  on  its  speed,  normally  significant, 

is  here  cancelled  by  the  effect  of  flow  divergence. 

“2  2 

When  P.  is  f(s  c)  the  curvature  is  0(e  ),  as  for  the  disturbed 
clans  flame  described  by  equation  (10).  We  would  therefore  expect  such 
large  flames,  if  disturbed  to  the  same  extent ,  to  be  described  by  a  modifi¬ 
cation  of  that  equation;  the  additional  terms  will  be  due  to  flow  divergence 
ar.i  undisturbed  curvature.  The  modified  equation  can  be  derived  by  formal 
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expansion,  as  was  the  original;  but  such  an  exercise,  although  reassuring, 
is  hardly  illuminating.  We  shall  instead  give  a  heuristic  derivation  that 
emphasizes  the  physics. 

The  general  effect  of  large-scale  (and  therefore  weak)  stretch  on  flame 
speed  was  identified  in  section  i*.5;  we  may  write 

w  =  i  -d+I)K+. .. ,  (18) 


where  K  is  the  Karlovitz  stretch  (^.3).  This  effect  is  the  origin  of 

the  second  term  on  the  right  side  of  equation  (6);  the  first,  corresponding 

to  what  is  normally  a  small  correction  to  the  result  (18),  must  be  retained 

when  i  is  close  to  -1.  In  the  present  context,  where  flow  divergence 

-1  2 

generates  a  stretch  H  of  order  e  in  the  undisturbed  plane  flame, 
equation  (6)  has  to  be  replaced  by 


>  _  -1 
W,  =  4F  +  so  +  R 

l  yy w  yy 


(19) 


I’ote  that  ?  still  represents  disturbance  of  a  plane  flame,  so  that  the 
description  is  valid  only  up  to  0(s  ‘ )  values  of  y. 

The  kinematic  expression  (3)  for  the  wave  speed  is  also  modified, 
because  the  flame  sheet  is  moving  in  a  non-uniform  velocity  field.  We 
find 

W.  =  -~.-y?r/H-?/R-y2/eR2-j€Ff,  (20) 

u  J  J 

to  sufficient  accuracy,  so  that  combination  with  the  result  (19)  now  yields 
+  -R , _ +  e?  +  r/P  +  y?  /?.  =  -y2/eR2  -  1/P.  (21) 

-  -  «*  J  V1  v  .■  J  J  J 


This  has  the  stationary  solution 


=  -y  /2eR 


(22) 


corresponding  to  the  undisturbed  circular  flame.  Replacing  F  by  F+FQ, 
so  that  F  nov  represents  disturbance  of  the  circular  flame  sheet,  and 
using  the  scaled  variables  (9)  yield 


F  +  +  UF  +  F 

t  n  nnnn  nn 


(23) 


Comparison  with  equation  (10)  shows  that  the  only  new  term  is  yF. 

> 

The  linearized  form  of  equation  (23),  with  F  set  proportional  to 
exp(oct+iky) ,  was  considered  in  section  5-^-  Figure  5-5  shows  that  curva¬ 
ture  is  a  stabilizing  influence,  but  that  instability  occurs  for 


Y  <  Yc  =  1/16, 


(2^ 


corresponding  to  a  supercritical  bifurcation  with  wavenumber 


k  =  /e/6 . 
c 


(25) 


~o  show  this  we  write 


Y  =  Yc-62,  F  =  of,  T  =  t/62 


(2b) 


and  expand  f  is  a  power  series  in  6.  In  the  usual  way,  we  find  that  the 
leading  term  is  of  the  form 


_  ik  y _ -ik  y 

fQ  =  A(t  )e  c  +  A(t  )e  c  , 


(27) 


;nen 


dA/dt  =  A-AA2/36 


(28) 


if  there  is  to  be  no  secular  term  in  the  second  perturbation  of  f.  The 
equation  describes  the  evolution  in  (slow)  time  t  of  the  amplitude 
initial  value  to  the  final  value  |a[  =  6. 

! 


some 
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An  examination  of  the  first  perturbation  of  f  reveals  that  the 
crests  of  the  final  stationary  solution,  as  viewed  from  the  burnt  gas, 
are  sharper  than  the  troughs.  Moreover,  the  flame  temperature  is 
diminished  at  the  crests  (as  for  the  cellular  flames  discussed  in  section 
5.3),  so  that  they  are  darker  than  the  rest  of  the  flame  sheet.  Sharp, 
dark  crests  are  a  universal  feature  of  cellular  flames  as  observed. 

3.  Flames  'Tear  a  Stagnation  Point. 

Equation  (23)  is  only  one  of  a  class  of  evolution  equations  that 
describe  cellular  flames  in  a  variety  of  circumstances.  An  unusual 
example  corresponds  to  a  flame  located  in  weak  stagnation-point  flow 

v  =  S(-x,y)  with  3  =  0(e2).  (29) 

Two  changes  have  been  made  in  the  notation  of  section  k.5:  the  x,y-axes 
have  been  rotated,  so  that  the  wall  is  now  x  =  0  and  the  ur.disburbed 
flame  at  x  =  x#  >  0  (cf.  figure  k.3),  to  conform  with  notation  already 
established  in  this  lecture;  and  the  strain  rate  is  now  8,  since  the  e 
used  there  has  been  conscripted  as  a  small  parameter  here.  As  for  the 
cylindrical  frame,  we  shall  be  content  with  a  heuristic  derivation. 

The  undisturbed  flame  experiences  a  stretch  6  so  that,  if  it  is 
displaced  by  an  amount  sF,  its  speed  is 


W,  =  Uf  +  eF  +  S 

1  yyyy  yy 


(30) 


c  sufficient  accuracy  (cf.  equation  (19))-  The  kinematic  expression  for 
;e  wave  speed,  corresponding  to  the  result  (20),  is 


„2 


■A  =  (Sx„-l)/e-r  -8  yFv  -  8F  -  ieBx.F? 

J-  *  t  J  j 


(31) 


'rusting  these  two  expressions  for  t-.'  gives  a  formula  for  x#,  namely 


•  f 


-6.9- 


2  2 
x*  =  1/e  y  +  e  with  y  =  8/e  , 


(32) 


and  an  evolution  equation  for  F: 


F  +  If"  +  kF  +  F  +  Y(nF)  =  0. 
■  *  n  nnnn  nn  n 


(33) 


Comparison  with  equation  (10)  shows  that  the  only  new  term  is  ytnF)  . 

n 

The  generalization 


F  +  J(V?r  +  W4?  +  V2F  +  y(nF)  =  0 

t  -  n 


(34) 


accounts  for  disturbances  that  vary  in  the  z-direction  also.  If  we  now 
consider  disturbances  independent  of  n,  this  equation  reduces  to  the 
earlier  cr.e  (23)  with  n  replaced  by  x,  =  /ez.  Setting  F  proportional 
to  exp (at+ikz )  ana  lir.eariz'ir^  therefore  leads  to  the  dispersion 
relation  (5.49)  and  hence  to  figure  5o-  The  bifurcation  analysis  starting 
with  the  transformation  (2c)  is  applicable,  so  that  for  values  of  y 
slightly  smaller  than  l/lo  there  will  be  a  stationary  structure  charc- 
terized  by  dark  ridges  pointing  towards  the  burnt  gas. 

This  phenomenon  has  apparently  been  known  for  many  years.  For  upward 
propagation  through  sufficiently  lean  hydrogen— air  mixtures  in  a  standard 
flammability  tube,  the  flame  cap  is  divided  into  a  number  of  bright  strips 
or  ribbons,  separated  by  dark  lines  (figure  4);  it  seems  probable  that  this 
is  the  axisymmetrie  analog  of  the  nominally  plane  flame  considered  here. 

The  straining  flow  is  generated  by  the  gravity-induced  convection  of  the 
light  burnt  gas  behind  the  flame  (see,  for  example,  Buckmaster  &  Mikoalitis 

A  iifferent  tyre  of  disturbance  (which  can  be  combined  with  the  previous 
one  without,  however,  adding  to  the  discussion)  corresponds  to 
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Substitution  into  the  linearized  version  of  the  evolution  equation  (33) 
gives 

dA/dx  =  A(k2e-2YT-UkUe-UYT-Y),  (36) 

vhich  has  the  solution 

A  =  AQexp[(l-e_2'T)k2/2Y  +  (e’^-l^-YXj  .  (37) 

Figure  5  shovs  that  any  disturbance  eventually  has  a  decreasing  amplitude, 
although  for  a  time  the  amplitude  increases  if  y  is  less  than  l/l6. 

In  the  limit  x  -*■  the  solution  (37)  tends  to  zero.  We  conclude  that 
the  flame  is  stable  to  this  type  of  disturbance. 

The  chaotic  cellular  instability  found  experimentally  for  weak  strain¬ 
ing  suggests  that  all  disturbances  should  grow.  Moreover,  in  the  limit 
Y  ■*  0  the  theoretical  results  in  lecture  5  predict  instability  for  all 
(small)  wavenumbers.  These  facts  are  at  variance  with  the  conclusion 
above,  which  prompted  Sivashinsky,  Law  &  Joulin  (1982)  to  provide  the 
following  explanation. 

If  the  nonlinear  term  is  retained  in  equation  (33),  harmonics  are 
continually  generate i  and  these  may  grow  during  part  of  their  lifetime, 
crovidir.g  a  mechanism  for  sustaining  the  overall  growth  of  the  distur¬ 
bance  into  instability.  Tumerical  computations  confirm  this  notion  and 
suggests  that  the  necessary  values  of  y  are  significantly  smaller  than 
1/1 C.  There  may  be  implications  for  the  hydrogen  flame  of  figure  k.  Away 
from  the  nose  the  rate  of  strain  will  be  diminished,  and  may  be  small 
enough  in  the  skirt  to  permit  instabilities  in  the  direction  of  the  flow. 
That  is,  the  ribbon  instability  cay  become  a  cellular  instability.  Inter¬ 
estingly  enough,  the  tails  of  the  ribbons  are  often  seen  to  break  up  into 
small  balls  of  flame  (figure  ^). 
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H.  Polyhedral  Flames. 

In  the  90  years  since  Snithells  first  observed  these  flames,  they  have 
become  established  as  a  familiar  laboratory  curiosity.  They  are  associated 
vith  tube  burners,  but  analogs  can  be  created  with  different  geometries: 

Mark stein  (196^,  p.  8l  )  designed  a  slot  burner  on  vhich  he  observed  a 
cellular  flame  behaving  essentially  like  an  unwrapped  (linear)  polyhedral 
flame.  In  particular,  the  cells  could  be  made  to  travel  rapidly  from  one 
end  of  the  slot  to  the  other,  just  as  the  polyhedral  flame  can  be  made  to 
spin.  Markstein's  photographs  of  the  flame  showed  that  the  traveling 
corrugations  are  saw-toothed  in  shape  (figure  6). 

It  is  the  propagation  that  distinguishes  polyhedral  flames  from  other 
types  of  cellular  instability,  so  that  will  be  the  focus  of  our  discussion. 
Since  the  left  stability  boundary  is  not  associated  in  any  obvious  way  with 
propagating  disturbances  (unlike  the  right  stability  boundary),  the  challenge 
is  to  uncover  a  mechanism  for  such  behavior. 

One  of  the  difficulties  with  polyhedral  flames  is  that  the  undisturbed 
flame  is  conical,  i.e.  non-piar.ar.  Following  Buckr.aster  (1933),  we  shall 
overcome  this  obstacle  by  adopting  a  nominally  planar  model.  Consider  the 
portion  .-13  of  the  burner  flame  that  is  located  near  the  rim  (figure  7a). 

The  flame  speed  varies  from  a  small  value  (perhaps  zero)  at  A,  to  a  value 
comparable  to  the  adiabatic  flame  speed  at  B.  This  portion  is  modeled 
by  a  plane  flame  with  some  intermediate  speed  and  standoff  distance  (figure  7b 
Perturbations  of  the  planar  configuration  are  permitted  in  the  y-direction, 
which  is  perpendicular  to  the  page  and  parallel  to  the  rim.  Corrugations 
that  arise  in  this  way  are  associated  with  the  corrugations  along  the 
entire  length  of  the  nominally  concial-shaped  flame  whose  base  is  being 
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In  the  context  of  the  weakly  nonlinear  theory,  perturbations  are 

governed  by  equation  (10),  provided  a  term  is  added  to  account  for  the 

presence  of  the  rim.  The  rim  is  a  heat  sink  that  anchors  the  flame  in  a 

simple  fashion:  an  increase  (decrease)  in  the  standoff  distance  reduces 

(increases)  the  heat  loss  to  the  rim,  thereby  increasing  (decreasing) 

the  flame  speed,  a  restorative  mechanism.  If  F  is  now  the  perturbation 

* 

of  the  standoff  distance,  this  effect  can  be  represented  by  the  modification 


W,  =  kF  -(l+2)F  +  q5 

i  yyyy  yy 


2 

of  the  equation  (6).  Here  q  is  a  positive  constant  of  order  a  ,  and 
the  factor  -(1+2.)  has  been  restored  (since  the  parameter  a,  which  was 
earlier  equated  to  it,  will  be  given  a  different  definition).  From  the 
kinematic  result  (8)  we  now  see  that  the  governing  equation  is 


v  +  Jf2  +  Uf 

t  y  yyyy 


(1+£)F  +  qF  =  0, 


(39) 


a  result  identical  to  the  curvature  equation  (23)  when  the  scaling  (9)  is 
undone  and  ?.  replaced  by  q. 

The  linearized  form  of  this  equation  was  considered  ir.  section  2,  but 
ve  shall  interpret  the  analysis  differently.  Rather  than  fixing  £  and 
determining  the  range  of  unstable  wavenumbers  for  each  q  (there  R  ^ ) , 
we  shall  fix  q  ar.d  determine  the  range  of  unstable  wavenumbers  for  each 
.  Thus,  with  7  proportional  to  exp(at+iky),  the  stability  boundary 
a  =  0  is  seen  to  be  the  curve 


K  _  p 

Uk  +  (1+Ok+q  =  0 


Tot  all  values  of  k  are  admissible,  however,  because  an  integer 
number  of  wavelengths  must  fit  around  the  burner  rim.  The  length  L  of 
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the  circumference  provides  the  definition  of  the  snail  parameter  e,  namely 


)  2/t2 
£  =  4tt  / L  , 


(ti) 


and  then  the  restriction  is 


k  = 


Vlth  ^  =  ,  K  =  1,2,3,... 


(1*2) 


The  positions  of  the  corresponding  points  on  the  neutral  curve  depend  on 

2 

the  value  of  q;  in  drawing  figure  8  ve  took  q  =  ll^c  .  The  corresponding 


values  of  2  then  lie  in  the  order  .  For  larger 

2  3  4  o  l  o  i 

values  of  q  there  are  more  points  on  the  lower  branch  of  the  curve:  in 

1  2 

general,  for  q  >  **M  s  there  are  M. 

Those  points  for  which  the  (reduced)  lewis  number  of  the  mixture  is 

less  than  correspond  to  unstable  modes.  Thus,  each  determines  a 

I?  N 

supercritical  bifurcation  corresponding  to  a  stationary  "-sided  polyhedral 

flame.  Such  unir.odal  bifurcations  are  essentially  the  same  as  those  consider 

earlier  in  the  context  of  flame  curvature. 

Both  1  and  q  car.  be  varied  in  ar.  experiment  by  changing  the  mixture 

strength  and  flow  rate.  Ve  have  seer,  that  changes  in  q  will  move  the 

points  corresponding  to  k..  along  the  neutral  curve,  altering  the  stability 

characteristics  of  the  flame.  Changes  in  £  will  do  so  too.  For  certain 

choices  of  q,  two  bifurcation  branches  merge,  i.e.  £,j  =  £^  for  some 

.  The  solution  or.  such  a  merged  branch  corresponds  to  a  spinning 

polyhedral  flame,  provided  =  211. 

The  most  satisfactory  case,  frcm  a  mathematical  point  of  view,  is 

2 

v  =  2,  "  =  1,  (i.e.  c  =  loe  ),  for  then  the  merged  branch  is  the  rightmost 
cr.e  and  presumably  is  accessicle  as  the  first  manifestation  of  instability. 
However  cr.e-  or  two-sided  polyhedra  do  not  fit  comfortably  on  a  circle,  so 
that  is  rot  a  physically  satisfying  choice.  The  objection  does  not  apply 


to  the  choice  !•!  =  6,  II  =  3  (i.e.  q  =  1296c  ),  but  then  the  two  branches 

originating  at  Z;  and  l..  lie  to  the  right  of  the  merged  branch,  and 
4  5 

our  elementary  analysis  can  provide  no  evidence  that  the  latter  is 
accessible.  Buchmaster  does  not  resolve  this  difficulty  but,  instead, 
argues  that  the  favorable  comparison  of  the  mersed-b ranch  solution  with 
the  physical  flame  is  good  evidence  of  accessibility. 

To  ensure  that  k.j  and  k^.j  give  the  same  Z ,  we  must  take 

q  =  l6>lV,  (2*3) 


and  then 


*H  =  *2R  =  -1"20rr  £' 


To  determine  the  solution  on  the  merged  branch,  we  perturb  q  and  £ 

o 

away  from  the  values  (43,2.1)  by  0(5"),  where  6  is  a  small  perturbation 


At  the  same  time  we  write 


F  =  5i,  o  =  t/c 


ar.d  expand  f  in  a  power  series  in  5.  The  leading  term  is  found  to 
he  of  the  form 

-  iX.y  _  _  -ik-.y  .  ik  y  _  .  -ik  y 

f-  =  A ( t ) e  “  +  A(t)e  ■'  +  3(t)e  **  +  B(t)e  “  , 


,  -  .  ?  2 
A/3t  =  -2kT-A3,  3B/3t  =  jkf.A 


if  there  is  to  he  r.o  secular  term  in  the  perturbation  of  f.  Partial 

derivatives  are  used  because  A  and  3  also  depend  on  the  slow  time 
•) 

ft;  evolution  :n  :his  scale  determines  the  ultimate  amplitude  of  the 


spinning  flame,  hut  we  shall  not  pursue  the  matter  here. 
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Equations  (1*7)  have  solutions  corresponding  to 
=  A„(/3  sin(wc+b±k..y)-cos  2(.J)t+$±kTTy)]  with  oj  =  k,“A  ;  (U8) 

U  \J  ii  14  L't  U 

here  A^  and  d  are  real  and  constant  on  the  t-scale.  These  are  waves 

traveling  in  the  negative/positive  y-direction  with  an  amplitude-dependent 

chase  soeed  k,TA„.  The  shape  of  the  wave  resembles  the  sawtooth  crofile  in 
u  *  , 

figure  6,  and  Bucksaster  has  argued  that  the  propagation  speed  is  consistent 
with  the  rapid  rotations  seen  in  experiments. 


5.  Other  Cellular  Flames. 

So  far  ve  have  been  concerned  with  the  evolution  of  the  linear 
instabilities  associated  with  values  of  2  slightly  less  than  -1. 

Various  additional  effects  were  incorporated  into  the  basic  nonlinear 
theory,  and  others  could  have  been  (Sivashinsky  1983).  Our  final  remarks 
are  concerned  with  values  of  2  slightly  greater  than  -1,  where  the 
linear  stability  of  the  flame  can  be  destroyed  by  hydrodynamic  effects. 

The  weakly  nonlinear  description  is  now 


+  Ut 


nnnn  *nn 


=°  F  (n  ,r  ) 

*  '(f  ‘ 


dn  =  0  with  y  = 


a-1 


n-n 


2te 


3/2 


(U9) 


Comparison  with  equation  (15)  reveals  that  the  sign  of  ?nn  has  been  changed, 
because  r.ow  the  definition 


e  =  1  +  2 


(50) 


is  needed  to  obtain  a  positive  parameter;  n  and  t  still  have  the 
definitions  (9b, c) .  Without  the  integral  term,  ?  ■*  0  as  t  -*■  °°  what¬ 
ever  the  initial  conditions  are,  corresponding  to  linear  stability,  as 
before,  the  hydrodynamic  effects  (represented  by  the  integral)  are 
ies~abilicir.g.  V.ieheison  ar.d  Sivashinsky ' s  computations  show  that  a 
pro-ressi/c  wave,  consisting  of  stationary  cells,  eventually  forms 
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provided  the  f lane  is  not  too  large .  For  large  flames ,  the  chaotic 

cellular  structure  first  found  in  section  1  reasserts  itself. 

The  shape  of  this  progressive  wave  satisfies  a  much  simpler  equation 

in  the  limit  y  ■*  0,  i.e.  for  significantly  larger  departures  of  l 

2 

from  -1  than  of  a  fro n  1.  Evolution  is  then  on  the  scales  yn,  y  x 

rather  than  n,T,  so  that  the  fourth  derivative  drops  out.  If  a 

nrogressive  wave  is  sought  by  setting  F  =  -V,  and  F  is  neglected 

t  nn 

(a  valid  step  where  the  curvature  is  not  large),  the  result  is 

an,  (51) 


o  00  ,  _  •  F  (n) 

jF  +  y^  (— —  -  —  )F  (n)dn  =  0  with  V  =  yf  J- 


n-n 


a  nonlinear  integral  equation  for  the  slope  F  . 

The  rough  solution  found  by  Sivashinsky  was  not  very  satisfactory 
ar.d  so  McCor.naughev ,  Ludford  &  Sivashinsky  (1983)  recently  integrated 
the  equation  more  accurately.  A  continuous  periodic  solution  is  shown  in 
figure  9;  the  solution  for  any  other  period  can  be  obtained  from  it  by 
scaling  n  without  changing  V  (there  is  no  preferred  wavelength  in  the 
linear  theory).  At  the  cusps,  F^  has  a  logarithmic  singularity,  which 
makes  the  structure  of  the  combustion  field  quite  different  from  that  of 
a  Bunsen  flame  near  its  tip,  for  example.  (Of  course  the  singularity  will 
be  smoothed  out  by  the  neglected  term  7  ,  as  Michel son  and  Sivashinsky ' s 

commutations  show.) 


V  =  l.Uy£ 


obtained  by  McCor.r.aughey,  Ludford,  and  Sivashinsky.  For  an  expansion 
ratio  of  c  =  3,  this  leads  to  a  flame  speed  that  is  1.6  times  the  plane 
adidacatic  value,  a  result  in  surprisingly  good  agreement  with  measured 
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values  (vhich  range  between  1.5  and  2).  The  theory,  which  assumes  o  -1 
to  oe  small,  is  not  valid  for  such  large  expansion  ratios,  but  nevertheless 
makes  an  accurate  prediction.  This  type  of  success  is  one  of  the  hallmarks 
of  a  great  theory'. 
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Figure  Captions 


6.1  Chrysler  emblem. 

6.2  Five-sided  polyhedral  flame. 

6.3  Numerical  calculation  of  cellular  flames. 

(Courtesy  G.I.  Sivashinsky . ) 

6.4  Flame  in  a  standard  flammability  tube. 

2  -2yt 

6.5  Curve  in  k  e  ,  y-plane  determining  sign  of  right 
side  of  equation  (36). 

6.6  Analog  of  spinning  polyhedral  flame  for  slot  burner. 

6.7  (a)  Behavior  of  tube  flame  near  rim  of  burner. 

(b)  Plane  model  of  (a)  used  to  describe  polyhedral  flames. 

6.8  Linear  stability  regions  for  polyhedral  flames,  with 
admissible  values  of  k  . 

6.9  Stationary  wrinkling  of  an  otherwise  stable  plane  flame 
due  to  hydrodynamic  disturbances;  possible  outcome  of 
Darrieus-Landau  instability. 
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